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Abstract 

The precise knowledge of the onset of the cross section a(e + e~ — * t + t~) at the 
threshold is necessary for improving the accuracy of determination of the r mass from 
the threshold measurements. The QED radiative corrections of relative order a and a 2 , 
additional to the well known Coulomb factor, are considered in the threshold region. 
The correction terms of order a 2 are calculated, which contain coefficients enhanced 
by large parameters. As a result it is argued that the known 0(a) corrections provide 
the accuracy of the description of the cross section close to 10 -4 , rather than 10 -3 
as claimed in a recent literature. Also analytical expressions are provided for some 
limiting cases of the corrections, previously calculated numerically. 



1 Introduction 



The production in the electron-positron annihilation of slow r + r~ pairs near the threshold 
provides a valuable tool for a precision measurement of the r lepton massQ. Such measure- 
ment is aided by that the threshold onset of the cross section starts with a finite step due 
to the Coulomb attraction between the produced r leptons||. Namely, when expressed in 
terms of the ratio R = a(e + e~ — ► T + r~)/a(e + e~ — > the familiar 'bare' threshold 

behavior 

Ro = v^f (1) 

in terms of the velocity v of each of the produced r leptons in the cm. frame, is multiplied 
by the Coulomb factor[] 

F. = i m / V n . (2) 
1 — exp(-ira/v) 

The product R F c clearly has a finite limit equal to 37ra/2 at v — > 0. The factor F c sums up 
all the graphs with exchange between the r leptons of Coulomb quanta and expands in power 
series in the parameter a/v, rather than in powers of the QED coupling a. For this reason 
at the near threshold energies, where v is not parametrically larger than a, the Coulomb 
interaction has to be taken into account exactly. The dependence of F c on a/v also implies 
that the higher order corrections can arise both from the QED radiative effects as powers 
of a as well as from the relativistic expansion through terms with extra powers of v 2 . The 
latter terms, being modified by the Coulomb effects, are generally equivalent to the ones 
with extra powers of a 2 . The leading QED radiative corrections of order a to the discussed 
process were analysed in Ref. [Q], where it was shown that these corrections arise from two 
sources: a velocity independent factor due to the form factor of the r electromagnetic vertex 
at the threshold, and the Uehling-Serber modification of the Coulomb potential due to the 
electron vacuum polarization, whose effect has a quite nontrivial dependence on v. It should 
be emphasized that the QED radiative corrections discussed here are those which are on top 
of the Coulomb factor F c . In other words all the terms of the form (a/v) n are accounted 
for exactly, and the discussed corrections are the terms of the form a(a/v) n : the 0(a) 
corrections, and a 2 (a/v ) n : the 0(a 2 ) corrections. 

Recently an attempt has been made|| at evaluating the 0(a 2 ) corrections, which could 
potentially contain large factors that would make them numerically significant. It was con- 
cluded in that the coefficients in front of the a 2 terms are not extraordinarily large and 
1 F r is also known as the Sommcrfcld-Sakharov factor. 
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that the 0(a) radiative corrections || are sufficient for describing the excitation curve for 
the t + t~ pair near the threshold with the relative precision not worse than 10 -3 . The cal- 
culation in Ref. || was based on essentially an adaptation of the NRQCD methods used 
to describe the QCD effects near a heavy quark flavor threshold (see the references in pj). 
The known results for QCD radiative effects within these methods are based on the MS 
renormalization scheme, and were used in such form in Ref. Q. It should be noted however 
that in QED, unlike in QCD, the renormalization is constrained by the simple requirement 
that the asymptotic behavior of the Coulomb potential at long distances is given by a/r 
with no corrections in terms of the physical fine structure constant a, tabulated^ as 
a^ 1 = 137.036. This property, inherent in the "on shell" scheme, is lost in the MS scheme, 
in which the Coulomb potential does receive radiative corrections at asymptotically long 
distances. Clearly this effect is entirely spurious, since it goes away once the QED constant 
in the MS scheme, otjjg, is expressed in terms of the physical constant a. However the 
necessity of including and subsequently eliminating the spurious terms somewhat compli- 
cates, if not obscures, intermediate calculations in the MS scheme. As an illustration of 
a difficulty arising in the latter scheme even at the level of the leading corrections it can 
be noticed that no distinction has been made in Ref. Q between the electron and the muon 
vacuum polarization in the Uehling-Serber type correction to the Coulomb potential. In 
reality, however, these two effects result in corrections to the cross section of completely 
different magnitude: the effect of the electron loop at the threshold is an enhanced by a 
factor containing ln(m T a/m e ) correction of order a, while that of the muon loop can rather 
be classified as an enhanced by the factor m T /m^ correction of order a 2 (or, alternatively, 
as an 0(a) correction, suppressed by m T a/m^). Numerically, at the threshold the former 
correction is almost 40 times larger than the latter 0. 

The purpose of the present paper is to demonstrate that the calculation of the 0(a 2 ) 
radiative corrections is significantly more transparent and simple in the "on shell" renor- 
malization scheme and to present an evaluation of the terms in these corrections, enhanced 
by parametrically 'large' coefficients. The latter 'large' coefficients include all the factors 
singular in the limit m e — > as ln 2 m e and lnm e , the already mentioned correction due 
the muon vacuum polarization, and also an enhanced 0(a 2 ) effect of relativistic corrections, 
which is proportional to the factor Inv (becoming In a at v <C a). As will be shown, the 
combined effect of the enhanced 0(a 2 ) corrections is only about 10~ 4 in terms of the rela- 
tive magnitude of their contribution to the cross section, while the non-enhanced terms are 
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proportional to (a/V) 2 rs 5 • 10 -6 , thus allowing one to argue that factually the corrections 
calculated in Ref. || already provide the theoretical accuracy of 10 -4 in the cross section, 
rather than 10 -3 as estimated in Ref. Q]. In practical terms, this implies that the theoretical 
accuracy is sufficient for a measurement of the r mass down to at least 0(1 keV), provided 
that similar experimental accuracy can be achieved in measurements at the r threshold. 

In Ref. II the effect of the 0(a) correction due to the vacuum polarization loop, both the 
electron and the muon, was presented in a form of a two-dimensional integral, which then 
was calculated numerically. This effect is also revisited in the present paper, and explicit 
analytical expressions will be given for the electron loop contribution in the limits v ^> a 
(but still !)<1), and v — ► (v <C a), as well as for the muon loop contribution, applicable 
for all values of the velocity below approximately m^/m T . 



2 Types of radiative corrections to a(e + e — > t + t ). 

Generally the QED radiative corrections in the actual cross section of the process e + e _ — ► 
t + t~ arise from the following sources^, 

i - radiation from the initial electron and positron, 

ii - vacuum polarization in the time-like photon, 

in - corrections to the spectral density p(q 2 ) = — | J2x(Q\jfi(~ °f the e ^ ec " 
tromagnetic current j M = (f 7 At r) of the tau leptons. 

iv - interference between the effects i - in which starts from the (relative) order a 2 . 

The actual cross section at the 'nominal' energy W = \fs in cm. of the electron-positron 
collision can thus be written in the form: 

rW 

a(W) = / r(W,w) |1 — U(w)\~ 2 a(w) dw + (interference terms) . (3) 

The weight function r(W, w) describes the radiation from the initial state@ and |1 — 11(^)1 _2 
is the factor for the vacuum polarizationfl in the time-like photon. These two effects are 
standard and are automatically accounted for in the data analyses, while the dynamics of 
the final state is encoded in the cross section a(w) = 87r 2 a 2 p(w 2 )/w 4 . The last term in 
eq.(|3]) arises from the graphs, where the lines of the initial electron and positron and of the 
r leptons are connected by more than one photon propagators. 

The subject of primary interest in the previous studies as well as in the present one of 
the discussed process is the spectral density p(q 2 ). However before proceeding to a detailed 
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discussion of the QED corrections in p few remarks are in order concerning the leading 
contribution of the interference. Clearly, this contribution can arise at the order a 2 (relative 
to the 'bare Coulomb' cross section), and is thus within the scope of the discussion in the 
present paper. 

One potentially possible contribution in that order could arise from the square of box- 
type graphs, where the r pair is produced through two photons. However two photons 
produce the r pair in a C-even state. For non relativistic heavy leptons the states can 
be classified in the standard (total spin) - (angular momentum) terms: 2S+l L with the 
production amplitude behaving near the threshold as v L (modulo the Coulomb effects that 
'convert' powers of v into powers of a at v ~ a). The production of the C-even S-wave 
state 1 S by the e + e~ is suppressed for chirality reasons by the factor m e in the amplitude, 
which for all practical purposes makes it totally negligible. The amplitude of production of 
the allowed by chirality C-even P-wave states 3 Pi and 3 p2 contains an extra power of the 
velocity, and thus the contribution of the box graphs to the cross section near the threshold 
is additionally suppressed by the factor v 2 |J . It can be noticed however that the amplitude 
for production of the 3 P2 state actually contains a somewhat enhancing factor lnt> M. Thus 
the relative magnitude of the correction due to the box graphs can in fact be estimated as 
5a /a ~ a 2 v 2 In 2 v ~ a 4 In 2 a, which is still far too small. 

The only other contribution of the relative order a 2 from the interference graphs can arise 
from the interference of the three photon production amplitude with the 'bare' amplitude 
mediated by one photon. The graphical representation of this contribution is shown as a 
unitary cut in the graph of Fig. la. We are interested here in the terms of the lowest order in 
the velocity, thus in the production amplitudes the velocity can be set to zero. One can verify 
that no singularity arises from the electron propagators in the limit m e — ► and thus that 
this term does not contain enhancing factors singular at m e — > 0. In other words, it does not 
contain factors with powers of ln(m T /m e ). This nonsingular behavior can be seen e.g. based 
on the Kinoshita-Lee-Nauenberg (KLN) theorem |T0|, |TT|j . Indeed, assume temporarily that 
the electric charge of the electron, Q e , and of the r lepton, Q T , are independent parameters. 
According to the KLN theorem there should be no infrared singular terms in the sum of 
the probabilities including the emission and absorbtion of soft photons, and this behavior 
is valid at arbitrary values of Q e and Q T . Expanding this total probability in powers of Q e 
and Q T one finds that in the order Q\ Q\ the only contribution to this probability in the 
limit v — > comes from the unitary cuts of the types shown in Fig.l, since at v — > there 
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are no cuts that would go across the photon line as well as across the r pair, and also the 
absorbtion of soft photons by the r leptons is irrelevant, since at v — > the heavy leptons 
neither radiate nor absorb photons. On the other hand the cut across the three photons 
as shown in Fig. lb does not contain infrared singularity in the limit m e — > 0: there is no 
singularity in the integrals over the energies of the photons, since the amplitude t + t~ — > 37 
has no such singularity in the photon energies^, while the collinear singularities also do not 
appear in the limit m e — > 0, since there is only one electron propagator factor per each 
photon. 




Figure 1: A representative graph (a) for the interference correction in the process e + e~ — > 
t + t~ . The sum of the unitary cuts in a and b contains no terms singular in m e at m e — > 
at the threshold. (Dashed lines show the unitary cuts.) 

The absence of photon radiation by the r leptons in the v — > limit also guarantees 
that no terms containing Inv arise from the graphs of the type shown in Fig. la. Thus this 
contribution can only appear as an a 2 correction in the cross section with just a numerical 
coefficient and containing no parametrical enhancement. As such this contribution should 
be taken into account in a complete calculation of the a 2 terms, which however is beyond 
the intended accuracy. 

Summarizing the previous discussion, one concludes that all the QED radiative correc- 
tions with relative magnitude of the first order in a and the parametrically enhanced ones 
of the second order, a 2 , are contained in the corrections to the r pair spectral density p(w 2 ), 
i.e. they originate from the source labeled as Hi above. 
2 This of course is known since the work |T5|. 
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3 First order radiative corrections. Effect of the elec- 
tron vacuum polarization 

In the lowest, 'zeroth', order the effective cross section a(w) can be expressed in terms of 
the nonrelativistic Green's function G(x, y, E) of the motion in the center of mass of the r 
lepton pair at energy E = w — 2mQ 

a(w) = ^-Im G(0, 0; m T v 2 ) . (4) 

In the absence of QED radiative effects the interaction between the r leptons is the Coulomb 
attraction, thus the Green's function is the well known one for the Coulomb potential V(r) = 
—a/r: G c (x, y, E). The imaginary part of the latter at x = y = is related to that of the 
free-motion Green's function, 

p 2 m exp(ip|x-y|) 

G (x,y;— ) = - , , , 5 

m An |x — y| 

by the Coulomb factor (Q): ImG c (0, 0; m r v 2 ) = F c Im G (0, 0; m T v 2 ), so that the leading 
order expression for a reads as 

tx 2 o? 1 

a ° = 1 1 T\ ■ ^ 

Lm* 1 - exp(-7ra/v) 

At the next level of approximation, i.e. in the first order in a, the corrections to o arise 
from two sources: from the so-called hard correction due to a finite radiative effect in the 
t electromagnetic vertex at the threshold, and from the modification of the Coulomb inter- 
action due to running of the coupling a, which is described by the Uehling-Serber radiative 
correction to the potential. The behavior of the Green's function at small separations x and 
y is determined by dynamics of the r leptons at characteristic distances r c ~ l/p c , where 
p c ~ m T v for v ^> a and p c ~ m T a for v ~ a and vCa. The hard correction to the vertex 
comes from distances of order l/m T , which are thus point-like on the scale of characteristic 
distances in the Green's function. Thus these two effects can be separated in terms of eq. @ 
[0,1 



as 



a{w) = (l- — ) Im G c (0,0;m T v 2 ) + 5 {1) G (0,0; m T v 2 ) 



mi V 7r 



<xo 11-—) (l + ^h(v)) . (7) 



71 ) V 3tt 



3 For a discussion see Ref. 101 and references therein 
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where the hard correction factor 1 — Aa/ir is well known in QED (see e.g. in the book ||13||), 
and S^G is the first-order correction to the Green's function due to the Uehling-Serber 
correction 5^V(r) to the Coulomb potential 

<jWG(x, y; m T v 2 ) = - J G c (x, r; m T v 2 ) 5^V(r) G c (r, y; m T v 2 ) d 3 

with 5^V(r) given by (see e.g. in the textbook |b 



j3 r ? 



The correction term h(v) in eq. (g) due to ImS^G can be found by considering the 
modification of wave function at the origin, 

5«V(0) = - J G c (0, r; m T v 2 ) 6V(r) tp c (r) d 3 r , (10) 

where ip c (r) is the S-wave wave function at energy m T v 2 in the Coulomb field —a/r: 

^ c (r) = C e~ ipr + i\2,2ipr) (11) 

with p = m T v, A = m T a/(2p) = a/(2v) and C = V'c(O)- Using the representation of the 
Coulomb Green's function in the form 

G c (0,r;p 2 /m T ) = -<^e^ jH e 2 ^ (i±^ dt , (12) 

and after the integration over r in eq. (fLOf) the result |J for h(v) at arbitrary v is expressed 
in terms of a double integral: 

«A , • „ „,_1\»A-1 



l + t\ iX (t + izxv- 1 ) ( 1 W 



/i = -2A Im / dt / dx — — i— ! — 7 ttt 1 + — ^ - — = — (13) 

Jo Ji V t J (t + l + izxv- 1 )^ 1 ^ 2xV x 2 1 ; 

with z = m e /m T . 

The integral in eq.flTBD can be readily calculated numerically 0] for an arbitrary relation 
between v and a. However it is instructive to have analytical expressions at least in the 
limiting cases: v ^> a and v -C a. In either case the discussed correction arises due 
to the running of the coupling a at distances r c ~ p" 1 , much shorter than the electron 
Compton wavelength. Thus the structure of term with h(v) can be readily understood by 
replacing the coupling a in the Coulomb factor F c by the effective coupling at momentum 
p c : a — * a (1 + 2 f- In J2iL ), thus finding 



= l-(l + 7ra/»>)exp(-7ra/r,) ^ ^ 
1 — exp(— 7r a/w) m e 
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In terms of this interpretation the subject of an actual calculation of h(v) reduces to finding 
p c in terms of m T a and m T v. The factor in front of \n(p c /m e ) in eq.(|14"D is obviously given 
by (a/F c ) (dF c /da). 

In each of the limiting cases the asymptotic behavior of p c is fixed up to constants: 
p c — > const ■ m T v at v 3> a, and p c — > const ■ m T a (with a different constant) at v — > 0. In 
the former limit of large v in order to find the constant one can use eq.(|8|) with the Coulomb 
Green's function replaced by the free one (eq.(|D), and also make use of the short distance 
limit of the Uehling-Serber correction: 

2q /, 1 5- 



*< Mr) = -^[l + ^(ln^-7*-|) 
r V Sir \ m e r bj 



(15) 



with 7e being the Euler's constant. The integral in eq.(|§) can then be readily done, and the 
result reduces to the replacement 

ln ^ ln 2m^_5 (16) 

m e m e 6 

in eq. fllil) at v ^> a (but still v <C 1). Taking also the u/a>l limit of the factor in front 
of the logarithm in eq. fllil) , one finds the expression for h(v) in this limit as 

i / m 7T / 2 Tfl T v 5 \ 

'*'->. = S (^-6)' (17) 



Compared to the result of the numerical integration in eq.(|T3D this expression provides the 
accuracy sufficient for description of the cross section with a relative error less than 10 -4 
down to v ~ 0.04. 

One can also find analytically the behavior of h(v) at v — *■ 0, using eq.(^) and the 
explicit form of the Coulomb Green's function at zero energy: 



ml a 



G c (0, r; E -> +0) = - ~ [F 1 (2 V ^^7) + i ^(2^^)] , (18) 

where J and K stand for the Bessel functions of respectively the first and the second kind. 
After inserting in eq.(|IT]) the correction to the potential in the form (^), the integral over r 
can be readily done, and the resulting expression for h(0) reads as 



Mo) 



m T a ( m T a\ (m T a 
exp K\ 



1 \ Vx 2 - 1 , 



m e x \ m e x / \m e x 
with K being the standard notation for the modified Bessel function of the second kind. 
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The integral in eq.(plj) can be evaluated using the presence of the large parameter 
m T a/m e rs 25.4. The integration over x can be split in two intervals: 1 < x < X and 
x > X, with X satisfying the conditions 1 < I < m T a/m e . After making the obvious 
approximations in these two intervals one finds the result as 



h(0) = In — !— + lE + - + O 

m P 6 



(20) 



which in terms of eq.flT^j) corresponds to the replacement 

P, 



In- 



m T a 1 
m VlE + ~ 



(21) 



Numerically, eq.(pO|) gives h(0) = 3.980 with a very high accuracy, given that the subsequent 
terms start with the fourth power of the small parameter m e /(m T a). 



4 Parametrically enhanced 0(a 2 ) vertex correction. 

The effect of the muon loop arising from its contribution to the Uehling-Serber correction to 
the potential can be evaluated by replacing m e — > m M in eq . (|i~3"D and considering there the 
limit z ^> 1. The resulting correction term h^(v) is then found to be essentially constant in 
v up to v w 2m A( /m r 0, and therefore in this region of v this term can in fact be approximated 
by /i( M )(0) from eq. (|I~9D . Using the latter expression and the fact that m T a/m^ <^ 1 one finds 



, . m T a f°° ( 1 \ \j x 2 — 1 , 9n m T a . . 

h ^ = ^l ( 1 + 2^) —^ dx= 32i^- < 22 > 

It can be noticed that the contribution of the muon vacuum polarization loop, as well 
as that due to the heavier states, to the discussed correction for the r pair spectral density 
at small velocity v of the r leptons behaves quite differently from the contribution of the 
electron loop. The reason for this difference is that a state contributing to the absorptive 
part of the vacuum polarization at q 2 = s modifies the interaction between the r leptons at 
distances shorter than ~ 1/ y^, thus the effect in the dynamics of the produced r leptons can 
be considered as local, provided that the velocity satisfies the condition v \fsjm T . Even 
for the muon threshold, \fs = 2m M , this condition is satisfied in all the region of interest for 
the velocity. The local effect does not depend on v in this region and can be calculated by 



4 This behavior also agrees with the numerical results of Ref. | 
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setting v — ► 0. Therefore in terms of the factorization formula (equations (^) and (0)) the 
effect of the muon loop and of higher states should rather be written as contribution to a 
correction, (a/7r) 2 A, in the vertex factor: 

. , 27r 2 a 2 / 4a a 2 A \ 0N . , 

aH = s- l + ImG0,fl;m T !; 2 . 23 

m% \ ix ix l ) 

Being interpreted in terms of the correction term A in eq. (p3j) the expression (^) corresponds 
to A( M ) = 37r 2 m T /(16m At ), which although being enhanced by the factor m T /m^, still gives a 
correction of only about 1.7 x 10 -4 in eq . (|23|) . 

The discussed correction due to the muon loop is only a part of the correction to the 
electromagnetic vertex of the r lepton at the threshold arising from the vacuum polarization. 
The latter correction can readily be found in full by noticing that the calculation can be 
formally reduced to a standard calculation of the one-loop vertex correction with a massive 
photon using the dispersion relation for the vacuum polarization P(k 2 ) with subtraction at 
k 2 = (the 'on shell' renormalization). Indeed, the photon propagator with one insertion of 
the vacuum polarization can be written as 

P ( fc2 ) - 1 f ds , (24 ) 

k 2 7T J s (k 2 — s) 

which clearly reduces the calculation to that with formally a massive photon with mass 

fJ. = Vs- 

At zero velocity, v — 0, the one loop vertex correction with a massive photon results in 
the following formula for the term A in eq.( 



with the function f(yfs/m T ) given by 



ImP( S ) JVAds (25) 



a \m 7 



(6 + 2^ + 2^- e 6 ) J—n arctan ^LJL + 2 ^ 5 ln^ - 3^ - 2^ 3 



(26) 



At y/s m T the function f(\fs/m T ) behaves as 



f(£) = ^-t + (£) . (27) 



Thus using the explicit expression for the muon vacuum polarization 



ImP(s) 1 ( A , 2m 2 \ 4 m 2 



a 3 



1 + \ 1 
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one finds from the integral in eq.(p5|) the singular in the ratio m r /m /J part of the muon loop 
contribution to A as 

37r 2 m T 8 , m T , , 

16 m M 3 m M 

The first term here clearly reproduces the correction given by eq.(^), while the second term 
describes the running a effect in the hard correction due to the fact that it comes from 
distances of order l/m T and thus should actually be written as — 4a(m T )/7i. 

The singular infrared behavior of the vertex correction makes it necessary to return to 
the discussion of the electron loop effect in the vacuum polarization contribution to the hard 
correction. For the electron loop, setting v — > in the calculation of the vertex correction is 
in fact not legitimate, since m T v (or m T a) is obviously not small as compared to m e , and the 
leading infrared term has to be accounted for within the Coulomb dynamics of the r lepton 
pair. This has been done in the previous section in terms of the correction factor h(v ). On 
the other hand the contribution to A resulting from the subsequent term of the expansion in 
eq . fl27|) and given by A( e ) = —(8/3) \n(m T /m e ) is perfectly correct, and describes the electron 
loop effect when a(m T ) in the hard correction is expressed in terms of the physical constant 
a. 

The contribution of the heavier than the /z 4 " [i~ pair states to the discussed correction A, 
in particular of the hadron vacuum polarization, thus contains no large logarithmic factors. 
Indeed, there is no 'logarithmic range' for such masses below the r mass, and the effect 
of the states heavier than r rapidly decreases with their mass (/(£) ~ 4 ln£ 2 /(3£ 2 ) at 
£ 2 = s/m T ^> 1). Thus keeping only the 'large' terms one can write the final expression for 
A as 

8 , m T 3n m T 8 , m T . . 

A = — - In — + — — - - In — 29 
3 m e lb m M 3 

Numerically these three terms almost cancel each other, resulting in an extremely small value 
of the vertex correction (a/V) 2 A w 1.0 x 10 -5 . 



5 Parametrically enhanced 0(a 2 ) corrections in the Green's 
function. 

In this section we consider the radiative effects of the second order in a in the Green's 
function in eq.(p3|). These effects arise from the second order iteration of the first order 
correction to the potential given by eq.(^), and from the second order radiative correction 
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5 {2) V{r) to the potential. It can be noted in connection with the latter effect that in the 
on-shell scheme there is no contribution to 8^ 2 'V(r) due to the diagrams with exchange of 
two photons between the r leptons, and all the corrections are given by the insertions in the 
single photon propagator, i.e. the two electron loop insertion and the correction in the one 
loop insertion. 

However, in order to find only those terms of order a 2 whose coefficients are enhanced 
by the second and the first power of ln(p c /m e ) there is no need for a detailed calculation of 
each these effects, and the result can be found using the KLN theorem. Indeed, the infrared 
singularity in ImG(0, 0, m T v 2 ) in the limit m e — > should be absent, provided that the result 
is expressed in terms of the effective coupling at the scale p c : a(p c ). Thus all the factors 
with ln(p c /m e ) arise through expressing a(p c ) in terms of the physical constant a. Up to 
the single logarithmic terms of the second order the latter expression is well known: 

/ \ 2a 2 , p c 4a 3 , 9 p c a 3 . p c 

« (Pc) = a + — ln^ + — In 2 ^ + — ln^ . (30) 
Sir m e 9tv z m e 2ir z m e 

Once p c is specified in terms of m T v and v/a from the results of the calculation in Section 
3, the discussed logarithmically enhanced terms can be found at any velocity in the region 
of interest. Writing the effect of these terms in the form of a multiplicative factor [1 + 
(a /it) 2 $(f )] in the cross section, one finds from the result in eq.(|20"D the expression for $ 
the limit v — > 0, where this correction is maximal, 

. - 4, 9 m T a /8 35 \ , m T a . . 

*(°» = 9 ln V + (9-» J+ M) ln ^7- (31) 

Numerically, the effect of this correction amounts to less than 0.5 x 10~ 4 . 

The function $(t>) decreases with the velocity, and at i; > a its behavior, as can be 
readily found from eq. (|16D , is given by 

, - . / n a\ /4, n 2m T v 13 , 2m T v\ . „. 

Hv) ^ ~ -In 2 1 In — . 32 

6 Relativistic corrections. 

Due to the presence of the Coulomb parameter, a/v, the relativistic effects, which formally 
arise from terms of order v 2 , get converted into corrections of order a 2 . Such corrections 
originate from relativistic corrections to the production vertex, and from corrections to 
the Green's function arising from the relativistic terms in the Breit-Fermi Hamiltonian. 
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A detailed study of these effects has been performed |T6|, |i~7j , |18f in QCD for the threshold 
production of heavy quarks in e + e~ annihilation. Since the results of this study in the 
order a 2 are not sensitive to the renormalization scheme, one can directly apply them to 
the discussed case of the r pair production near the threshold. The relativistic correction 
contains as 'large' parameter lnp c /m T , and keeping only the enhanced terms, it can be 
written as the following multiplicative factor in the cross section 



, 2 2 

1+ r 



1 fa 
In Re ^ i-i — 

v V 2v 



(33) 



where \I/ stands for the digamma function. 

Numerically, this correction term reaches its largest value [2 a 2 ln(2/a)]/3 ~ 2.0 x 1CT 4 
at v — and slowly decreases with velocity (e.g. the numerical value of the correction at 
v = 0.1 is 1.7 x 10~ 4 ). 

It can be also noticed in relation to the expression in eq.(^) that being multiplied by the 
Coulomb factor (0) it produces a double pole in the cross section at t> = in a/2 corresponding 
to the n-th bound 3 Si level of the r + r~ pair. This double pole correctly reproduces the Breit- 
Fermi relativistic correction to the energy of the bound state jnj. On the other hand the 
singularity of the digamma function at v — » is canceled by the logarithmic term. These 
two observations are in principle sufficient to reproduce the structure of the discussed part 
of the relativistic correction and the coefficient in front of it. 



7 Summary and concluding remarks 



We are now ready to collect all the discussed terms into one formula and to write down 
the complete expression for <j(e + e - ) — > t + t~ near the threshold including all corrections of 
order a and all the parametrically enhanced terms of order a 2 : 



a 



7r 2 a 3 



2m 2 . 1 
2tt 2 



exp(— na/v) 



4a 2a , , . 
1 + — h(v) 

TV 37T 



$(v) - 



m T 3ir m T 8 m T 
In 1 — — In 



16 m. 



m. 



1 fa 
In Re * [-1 — 

v V 2v 



3%) 



(34) 



where the function h(v) is given by the equations (0), QITD, and ( pD|) , and $(w) is defined 
in the Section 5 (cf. the equations (PT|) and (§^)). The last term with h{y) in the 0(a 2 ) 
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correction in eq.fl34]) obviously arises from the interference of the first order vertex correction 
with the electron loop effect in the first order (cf. eq.(|7|)). 

The relative magnitude of the discussed QED radiative corrections decreases with v. Thus 
their combined effect can be majorated by their value at v — 0. Writing the numerical values 
at v = of the individual terms in the inner curly braces in eq.(]33j), the 0(a 2 ) term can be 
evaluated as (a/vr) 2 (8.41-21.75 + 31.12-7.52 + 36.93-10.61) = 36.58 (a/vr) 2 « 2.0 x 10~ 4 . 

It can be also noticed that the 0(v 2 ) relativistic term, present in eq.([l]) was intentionally 
omitted throughout the previous discussion as well as the total energy dependence of the 
cross section arising from the timelike photon propagator producing the r+r~ pair. This 
is justified, as long as non enhanced 0(a 2 ) corrections are also ignored, in the range of 
the velocity where v and a are considered as being parametrically of the same order, and 
where the discussed effects of the interaction between the r leptons are most interesting. 
For practical reasons however it is desirable to have a description that interpolates with a 
high accuracy the cross section between the 'Coulomb' region of small v and the relativistic 
region. An interpolating formula of this type was suggested in Ref.[Q in terms of the well 
known full expression for the cross section up to order a (but without any summation of the 
Coulomb terms): 

(7o ( e + e ^ T + T ~) = ^! v ( 3 _ J) ^ + SL S (vf) , (35) 

where S(v) is given by 

S(V) = ~ { ^ + ^ + ln (¥) ln (S) + 2 Li2 fe) + 2 Li 2 {^) ~ 

2 Li 2 (1=2) - 4 Li 2 (v) + Li 2 (v 2 )} + v 2 ) -3v + ln (g) + (36) 

Qv ln (1±2) - Av Inv + \v { j0±) 

with Li 2 (x) = — J(fln(l — t)dt/t = Y^=i xn / n2 - The interpolation formula,which does 
include the summation of all the Coulomb terms as well as the 0(a) correction to them, has 
the form 

^(e+e- ^ r+r") = ^— i; (3 - t; 2 ) F c 1 + -S » - + tt" ■ (37) 

6S \ 7T 2f 07T / 

Near the threshold this formula correctly reproduces the Coulomb enhancement as well as 
the first radiative corrections to it. At all v the interpolation correctly reproduces the terms 
of the first order in the expansion in a. From these observations it was concluded that the 
higher corrections to eq.(|3"7]) are uniformly of order a 2 at any velocity with no enhancement 
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by inverse powers of v near the threshold. From the calculations in the present paper it is 
seen that the corrections of the relative magnitude 0(a 2 ) in fact receive a moderate logarith- 
mic enhancement near the threshold which however does not exceed 2 x 10 -4 numerically^. 
I believe that this accuracy should be quite sufficient for all foreseeable practical purposes. 
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